COMPLEX BODIES WITH MEMORY: LINEARIZED SETTING 



PAOLO MARIA MARIANO AND PAOLO PAOLETTI 



Abstract. The mechanics of complex bodies with memory effects is discussed 
in linearized setting. The attention is focused on the characterization of free 
energies in terms of minimum work and maximum recoverable work in the 
bulk and along a discontinuity surface endowed with its own surface energy, a 
surface internal to the body. To this aim, use is made of techniques proposed 
by Del Piero. Consequences of the Clausius-Duhem inequality are investigated 
for complex bodies with instantaneous linear elastic response. 



1. Introduction 

Some materials display sensibility to the past history of their present state: they 
are called the materials with memory. A paradigmatic (although special) example 
of description of hereditary behavior is the standard viscoelasticity, even in small 
strain regime, for example the one described by the rheological models by Newton, 
Maxwell or Kelvin. The pioneering work by Volterra [42l [43] opened the way to 
the analysis of complicated constitutive structures built up on the histories of state 
variables (basically their graphs in time). 

The topic received careful attention in the late 1950s, the subsequent decade 
and more. Researches were developed with the aim of establishing a theory of 
linear and non-linear behavior of materials with memory, especially in the case of 
fading memory [311 El IlHl [3 [6] . All these works generated specific developments 
and general effort toward ordering of the existing results in a clear and systematic 
way. The articles on this topic are manifold and it is difficult to list all of them. 
Essential examples are[2Il[71[a[ni[5Hl[31[Il[Il[ini[ni[IHl[ISl[ini[Il[lIl[ll 
[24[ [25l [26l [36l [38] , all dealing with the description of memory effects in Cauchy 
bodies that are those bodies the morphology of every material element of which is 
described by the sole place that its centre of mass occupies in space. 

Essential questions have been tackled in the representation of the mechanical 
behavior of the materials with memory. Some of them have foundational nature: 
(i) the meaning to be given to the notion of state, (ii) the definition of appropriate 
free energies, (in) the nature of the chain rule when functionals of histories are 
called upon, (iv) the correct use of Clausius-Duhem inequality to find a priori 
constitutive restrictions. 

In particular, in linear viscoelasticity items (i) and (ii) have been successfully 
tackled in [3T] and [25] : viscoelastic material elements of Cauchy bodies have been 
described as thermodynamic systems in the sense of [12] . 
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Here, techniques proposed by Del Picro in [21] are used extensively to address 
the mathematical modelling of the behavior of complex bodies with memory in the 
general model-building framework of the mechanics of complex bodies (different 
analyses on the specific case of micromorphic materials with memory are available 
in [^[30]). 

Bodies are called complex when changes in the molecular or crystalline texture 
at various microscopic scales (substructure) influence the macroscopic behavior 
through peculiar actions. Examples are manifold: liquid crystals, ferroelectrics, 
quasicrystals, nematic elastomers, magnetostrictive solids etc. Although all these 
examples are referred to a variety of phenomena, common essential features can be 
referred to a unique abstract model-building framework for the mechanics of com- 
plex bodies (see [21 IBS', "34]). Such a framework unifies in a single format the existing 
models of special classes of complex bodies and is a flexible tool for analyzing new 
materials. 

Some complex bodies, as for example relaxor ferroelectrics [H I44j , exhibit mem- 
ory of substructural events. Such a circumstance motivates the analysis of the 
mechanics and thermodynamics of bodies with memory. The setting is the linear 
one for the sake of simplicity. The representation of the substructural morphology 
of bodies is maintained abstract in order to include a number of special cases as 
large as possible. 

Free energies are characterized in terms of minimum work and maximum recov- 
erable work in the bulk and along a discontinuity surface endowed with its own 
energy. Consequences of the Clausius-Duhem inequality are investigated for com- 
plex bodies with instantaneous linear elastic response. 

2. Kinematics of complex bodies 

2.1. Generalities. A regular region B in the ambient space M'^ - regular in the 
sense that it is a 'fit region' or, more simply, an open set with Lipschitz boundary - 
is selected to host a body in its macroscopic reference configuration. Every point x 
from B represents a material element. It is assumed that subsequent configurations 
are achieved by means of differentiable bijections (transplacements) from S to a 
copy R'^ of the ambient space, obtained by means of an isomorphism i : R'^ ^ R'^. 
Maps X I — > y := y (x) G R"^, x G B, are then defined. Their spatial derivative is 
indicated by F := Dy (x) £ Horn {TxB,Ty(^)Ba) , with Ba ■= y{B), and is such 
that det > 0. Comparison between the metric g in the actual shape in Ba and 
the natural metric j in B allows one to measure crowding and shearing of material 
elements. The tensor E :— ^ {y'^g — 7), where y'^g is the pull-back of g through y 
given in components by {y'^g) = PjigijFs^ then a true measure of deformation: 
it vanishes under global rigid transplacements. Motions are then time parametrized 
families of transplacements: 

(a;,t)i — > y:=y{x,t) e'M?, x e B, t e [0,d]. 

Sufficient smoothness in time t is presumed. The macroscopic velocity is defined 
by y := -^y {x,t) in the referential description (that is as a field over the tube 
B X [0, d]). The standard kinematics of deformable bodies is represented this way. 
No geometrical information on the material texture at scales lower than the macro- 
scopic one (substructure) is commonly added in the description of the morphology 
of the body under scrutiny. 
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When materials display acute sensibility to phenomena at minute scales, a rep- 
resentation of such phenomena, combined with the description of the macroscopic 
behavior is required. The standard kinematics is then enriched. A descriptor of 
the material texture (a morphological descriptor of the substructure) is assigned to 
each point. It is selected in a set Ai. A field 

{x, t) I — > v := V {x,t) M, X e B, t e [0, d], 

is then introduced: differentiability in space and sufficient smoothness in time are 
presumed. The time rate of change of the morphological descriptor field in referen- 
tial description is ;> :~ {x,t) G T^(^x.t)M, the spatial derivative is indicated by 
N -.^ Diy{x) € Horn {T^B, T^(^)M). 

In selecting a specific morphological descriptor v of the material texture (a pa- 
rameter also called a fabric tensor or an order parameter) one choose the prominent 
geometrical features of the minute world inside the generic material element to be 
described, transferring their peculiarities at gross scale. However, the construction 
of the essential structures of the mechanics of complex bodies requires only that 
Jvi be considered as an abstract finite-dimensional differentiable manifold. Spe- 
cific geometrical property of M have often a clear physical meaning. A metric is 
the basic ingredient for the representation of the independent peculiar kinetic en- 
ergy pertaining to the substructure, if such energy can exist in special cases (see 
[3]). A connection allows one to represent contact substructural interaction (mi- 
crostresses) and to decompose in invariant way them from the self-actions occurring 
in each material element (see [39]). For these reasons, the specific nature of M. is 
left unspecified in the developments below. The subsequent results then hold for a 
wide class of complex bodies. 

2.2. A discontinuity surface. A surface 

E {x e clB, fix) = 0}, 

with / a smooth function, is selected in B. At x the normal m to E is defined by 

_ V/ (^) 
|V/(x)|' 

and orients E locally. Notice that by such a definition the normal is considered as 
a CO- vector. The projector over E is the second-rank tensor {I — m® m), with / 
the identity. Let x i — > a (x) be a differentiable field over B. Its surface gradient at 
x G E is given by V^a := Va [I — m® m). The trace of V^a defines the surface 
divergence of a at x, namely Div^a = trV^a. 

Assume that x <—>■ a (x) takes values in a linear space. If it is piecewise differen- 
tiable and suffers a bounded discontinuity over E, its jump [a] across E is defined 
by [a] := — at every a; G E, that is by the difference between the inner and the 
outer traces of a at E, defined by the limits a* := limj^o+ a {x i em). The average 
(a) of a across E at every x is defined by 2 (a) :— a'^ + . For every pair of fields 
fli and a2 with the same properties of a, the relation [aia2] — [ai] (02) + (ai) [a2] 
holds if the product 0102 is defined in distributive way. 

It is assumed here that both x 1 — > F and x 1 — > N are discontinuous across E, 
while the field x 1 — > ly is continuous there. The symbols F and N denote the surface 
gradients of deformation and the morphological descriptor, respectively. They are 
defined by 

F := (F) (/ - m (g) m) , N := (TV) (/ - m (g) m) . 
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2.3. Linearized kinematics. Convenience suggests the introduction of th.e dis- 
placement field 

(x, t) I — > u := u (x, t) = y {x, t) — i (x) , (x,t) £ B x [0,d] . 

The spatial derivative of x i — > u is indicated here hy W := Du{x,t). One gets 
obviously y (x, t) = ii (x, t) and F ~ I + W , with / the second rank identity tensor. 
The condition 

\W\ « 1 

defines the infinitesimal deformation regime. The deformation tensor E is then 
substituted by its linearized part e :— symW. Moreover, no distinction is also 
made between B and Ba in the sense that y « m at every x in B. No distinction 
is also made between x i — > v and y \ — > Va '■= v ° y^^ ■ Along the discontinuity 
surface S, the surface displacement gradient is indicated by 

W:= (VK) {l~m®m). 

The linearized kinematical setting justifies the mixed use in the same context of 
symbols adopted elsewhere for distinct actual and referential measures of interac- 
tions, as it is made in the ensuing section. 

3. Power and balance of interactions 

3.1. Classification of the actions. Relative changes of places between neigh- 
boring material elements generate standard actions represented in the Lagrangian 
description (that is as a field over B) by the first Piola-Kirchhoff stress P S 

Horn (t*B, Tyf^^^Bc^ and the vector of body forces h G R^*. Substructural events 

may occur within the material elements even when the material elements themselves 
are frozen in space. Inhomogeneous substructural changes in space generate new 

contact actions measured by the so-called microstress tensor S G Horn (t,*B, T*f^^-^A4 

External bulk fields can act directly over the substructures (magnetic and/or elec- 
tric fields, or some other radiative fields). They are represented by the co- vector 
P G PJ(j.)A^, at each x. All these interactions contribute to the expression of the 
power of all external actions on a generic part of the body, namely on any subset 
b of S with non vanishing volume and the same regularity properties of B itself. It 
is said that a generic part b crosses E - in this case it is indicated by bs - when 
dbT, n E is a simple closed curve where the normal n is defined as a vector from 
T*S at all X G 9bs H S where the normal to 9bs exists; in particular n at a given 
X belongs to the tangent plane to S at the same point. 

The surface S can be considered as a model of a material layer with vanishing 
thickness. In this case it is called a structured surface and it is assumed that it 
can carry standard and substructural surface actions, the former represented by a 
surface stress T, the latter by a surface microstress §. 

For a generic part bs crossing S, the explicit expression of the power is the given 
by (see [33]) 

Pb^^* {y,iy) : = / {b-y + P-u) dx+ ( [Pn ■ y + Sn ■ u) dU^ + 
+ [ (Tn-(7))+Sn ■(;>)) dn\ 
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A Lagrangian representation is used in the earlier formula. The link with the 
Eulerian (actual) description - a link given by the standard Piola transform - is 
recalled later. Then everything is reduced to the linearized setting. 

3.2. Observers. An observer is intrinsically a representation of all geometrical 
environments which are necessary to describe the morphology of a given body and 
its motion. 

The setting discussed here then incudes the assignment of atlantes over the 
reference place B, the ambient space M.^, the interval of time and the manifold of 
substructural shapes M. Changes in such atlantes are changes in observers. Amid 
them the interest is focused here on synchronous changes in observers - the ones 
leaving invariant the representation of the time scale - which evaluate the same 
reference place. In this sense only changes in W' and A4 are accounted for. 

The ambient space is altered by the action of the group of diffeomorphisms 
onto itself, namely the group Dif f{M.^ ,M.^). Its action has infinitesimal generator 
coinciding with the vector field which assigns to each point the vector ^ \s=o, 
where fg is a point selected over a smooth curve s i — > fs , s 6 in Diff{M.^, R"^) 
such that fo = identity. The parameter s can be identified with the time. 

Since the material substructures are in fact placed in space, changes of frames 
in M"^ alter in principle the geometry of the substructures and their consequent 
representation over A4. There is exception when v represents only a generic 
property of the material substructure not associated with its geometry in space. 
Besides this circumstance, one may presume the existence of an homomorphism 
h : Z)i//(K^,M'^) G, with G the Lie group of diffeomorphisms of M onto 
itself and h mapping the identity in Di//(R'^, R"^) to the identity in G. The 
curve s i — > is then generates a curve s i — > gs := ft-(fs) over G, and the cor- 
responding infinitesimal generator of the action of G over A4 is then defined by 

Cm (^) •= ^ \s=o = ^^ds^ ls=o (^^^ related discussions in [S^ and [35]). 

Changes in observers generated by the group SO (3) of the proper rotations, a 
subgroup of £)i//(R^, R'^) are specifically under scrutiny. For f\q an element of the 
Lie algebra so (3), g G R'^, one writes the corresponding {v), obtained through 
h, as the product A {v) q with A (v) e Hom{R^, Ti,M). 

By indicating by y* and v* the pull-back in the frame of the first observer of the 
rates evaluated by the second observer), one gets 

y* ^y + q^{y-yo) 

where yo is an arbitrarily fixed centre of rotation in the ambient space, and 

i>* ^ V + A{v) q. 

Here s is identified with the time. 

3.3. Invariance and its consequences. 

Axiom 1. At ( dynamic) equilibrium V^^^ (y, v) is invariant under rotational changes 
in observers. 

Theorem 1. (i) If for every bs the vector fields assigning the values an and A*Sn 
are defined over c^fas and are integrable there, the integral balances of actions on bs 

I b dx+ I Pn dU^ + I Tn dU^ = 0, 
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/ {{x-xo) Ab + A*P) dx+ {{x ~ xo) APn + A*Sn) dH'^+ 

+ [ iiy- yo) A Tn + ^*§n) dH^ = 0, 

(ii) Moreover, if the tensor fields x i — > P, S are of class {B\Yi) and are also 
continuous over the boundary of the body, then 

DivP + b = Q, 

and there exist a co-vector field x i — > z G T^.j^.-jA^ such that 
skw (PF*) = ie {A*z + {\7A*)S) 

and 

DivS - z + /? = 0, 

with z = Z1+Z2, Z2 G KerA*in the bulk. Additionally, if the tensor fields x 1 — >, T, S 
are of class (S) along the surface E and are also continuous along its boundary, 
one gets 

Div^T+[P]m = 0, 
and there exists a co-vector field x 1 — > 3 G T^(^^)M., with a; G S, such that 

skw (TF*) = ie {A*i + (Vs^*) §) 

and 

Div^S - 3 + [S]m = 0. 

(Hi) If the rate fields {x,t) 1 — > y{x,t) G K."^ and {x,t) 1 — > z>(a;,i) G T^(^x)^ 
differentiable in space, the local balances imply 

where 

VI"* (y» := ( [P-F + z-u + S-N) dx+ [ (T • F + 3 •!> + §• N) dH^ . 

"P™* {y, v) is called an inner (or internal) power, e indicates Ricci's alternating 
index. 

If one enforces Axiom 1 with a requirement of invariance with respect to the 
action of the semi-direct product k 5*0(3), rather that calling upon only the 
action of SO (3), a proof of Theorem 1 can be found in [33j . The weaker requirement 
here imposes a change in the proof. By using the Axiom 1, in fact, one first obtains 
only the integral balance of moments, then one first exploits the arbitrariness of 
the centre of rotation and substitutes j/o with j/q + w, with w an arbitrary vector 
depending only on time. If one subtracts the integral balance of moments from 
the resulting equation (the one obtained by the substitution y^ \ — > yo + w) and 
the arbitrariness of w also imply the integral balance of forces. Pointwise balances 
follow by the standard use of Gauss theorem (see also remarks in pHl[55] ). 

In the Eulerian representation the balance equations become 

diva + ba = 0, 
skw (cr) = -e [A* Za + [gradA*) So) , 

divSa - Za + /3a = 0: 
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in the bulk and 

div^Ta + [o-]ma = 0, 
skw (T„) - ie {A*ia + igradsA*) §a) , 

diVY:Sa - ia + [Sa]ma = 0, 

over E. Moreover, by indicating by v and v the values the values of the actual 
(Eulerian) representation of the (differentiable) velocity fields {y,t) i — > v {y,t) and 
(y, t) I — > V (y, t), one also get the equation 



{ba ■ V + P ■ v) dy + / (cTOa • V + Safla ■ v) dH^ + 

(TaHa ■ {v) + SaHa • (u)) dH^ = / (cT • ^mdu + Za • u + 5a • Qradv) dy+ 

Jy(h) 

+ / (Ta • gradsw +3a • w + Sa • gradsw) dH^ . 



ly(b)ny(T.) 

In previous formulas ria, ma and ria are the counterparts of n, m and n in Z?a. The 
differential operators dw, dws, grad and grad^. involve derivatives with respect 
to the coordinates in the current macroscopic placement. The index a means 
'actual'. Moreover, the actual measures of interactions are obtained by means of 
the standard Piola transform as 

6a := (detF)"^6, cr := (det F)"^ PF*, 

Za := (det Py^ z, (3^ (det /3, Sa (det P)~^ SP* , 
Ta := (detF)"^TF*, 3,, := (detF)"^ 3, Sa := (detF)"^ §F*. 

The proof of the Piola transform can be found on any textbook in nonlinear contin- 
uum mechanics. Less popular is its counterpart on surfaces embedded in a body: 
for the relevant proof see [29] . 

The Piola transform implies that, in infinitesimal deformation setting, referential 
and actual measures of interaction in the bulk and over the surface E coincide as 
\W\ and |W| tend to zero. 

The ensuing sections are just restricted to the infinitesimal deformation setting. 
Thus, by taking into account the substantial coincidence of referential and actual 
measures of interactions in the linearized setting (as remarked in earlier comments), 
although the Cauchy stress a is used, the index "a" in the other actual measures 
of interaction is omitted for the sake of conciseness. 

4. Linear constitutive structures 
Constitutive structures of the type 

a = a{W,v,N) , z^z{W,v,N), S ^S{W,v,N), 
in the bulk and 

T = T(W,z^,N), 3=3(W,i^,N), § = §(W,j^,N), 

on the surface E can be selected for elastic complex bodies. The entries of the 
previous constitutive structures are instantaneous value. Linearization of them 
about a pair (u,i>) requires the embedding of 7W in some linear space isomorphic to 
M'^ for some k. This embedding allows one to consider the space of pairs of maps 
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{y, v) as an infinite-dimensional manifold modelled over a Sobolev space: the use 
of the Frechet derivative in the linearization procedure then follows. 

It is then assumed that M. is endowed with a Riemannian metric and the 
relevant Levi-Civita parallel transport. Such structural assumption is constitutive 
in its essential nature. By Nash theorem an isometric embedding of in a linear 
space is then always available but it is neither unique nor rigid. The selection of 
an embedding plays the role of a constitutive ingredient of every special model. 

Under these conditions, constitutive equations expressed by linear operators L^'^, 
namely 

cr = L('") {W, iy,N), z = L(^) {W, u,N), S = L^'^^ {W, v, N) , 

in the bulk and 

T = L(^) (W, zy, N) , 3 = L(3) (W, u,N) , § = L^^' (W, z/, N) , 

on the discontinuity surface, make sense. 

The procedure discussed here holds also when the measures of interaction depend 
not only on instantaneous values of the state variables but also on their entire 
history. In this case memory effects can be accounted for. Viscosity come into play. 

5. Characterization of the bulk free energy in terms of work: 

VARIATIONS on A DEL PlERO'S THEME 

5.1. Histories. At a point x, a history is a BV right continuous map 

H:R+ ^ M3X3 X R'^ X Mk^s 

such that, for s e M"*", 

H{s) = {W{s),iy{s),N{s)). 

Its restriction over an interval [r, p) is called process and is defined by 

K^is) := {Wir + s),i'{r + s),N{r + s)) , < s < p - r. 

As shorthand notation, Kp indicates a process when it is of the type K^. The 
symbols F and 11 denote the spaces of histories and that of processes, respectively. 
Of course, here 11 C F. Processes prolong histories. Given H, the history 

is called prolongation of H by means of the process Kp. It is assumed also that 

Kp{p)- =H{0), 

where Kp {p)~ := VimsypKp (s), to assure differentiability in time, a property nec- 
essary for later use. 

Along S, a surface history 

s I — > H(s) := (W(s), v{s), N(s)), 

can be defined when the map x i — > y{x) is continuous across the surface. Of course, 

since Ad is embedded in a linear space, the average {u) makes now sense so that 
one may consider also ]HI(s) to be coincident with (W(s), {v) (s),N(s)). The results 
collected below hold also in this case. 
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5.2. History dependent measures of interaction and equivalence of his- 
tories. In the earlier notes, the state variables have been indicates just formally. 
When memory effects are accounted for, the notion of state requires careful defini- 
tion because different equivalence relations between pairs of histories exist so that 
the notion of state appears to be rather natural in terms of equivalence classes. 
Such a question has been tackled variously (see |22l [211 ES]) on the basis of the 
abstract approach to thermodynamics proposed in [12] (see also [13 1 l40 l l4T]). 

In what follows, the point of view developed by Del Piero in [3T] for linear simple 
bodies with memory is adapted to cover linear complex bodies displaying memory 
effects at macroscopic and microscopic scales. The aim is of (1) characterizing states 
and (2) deducing the main property of the free energy of such complex bodies. 
A concrete example of such bodies is the one of relaxor ferroelectrics. Butterfly 
loops in the diagrams of strain versus applied electric fields indicate the presence 
of memory effects [il] . 

Linear constitutive structures in the bulk are here assumed to be of the form 



a(F) = G^w{0)WiO) + G^AO)iy{0) + G^NiO)NiO) + 

+ ^ [CaW is) W (s) + (s) V (s) + G^A, (s) N (s)) 



S[H) = Gsw (0) (0) + Gs, (0) J. (0) + GsAT (0) (0) 



[Gsw (s) W (s) + Gsu {s) V (s) + GsN (s) N (s) 



z{H) = G,M/(0)W^(0) + G,,(0)i.(0) + G,Ar(0)A^(0) 







(g.w (s) W (s) + G,. (s) V (s) + G,N (s) N (s) 



ds. 



ds. 



ds. 



where the G^s's are the so-called relaxation functions, tensor functions (taking 
values in different tensor spaces) that are assumed to be Lebesgue integrable in 
time: they are absolutely continuous and the limit G^b(oo) := liixis^+oo Gab {s) 
exists. Of course, the indexes A and B run in {a, z, S} and {W, A^} respectively. 

Notice that the integrals above are well defined because M. is considered embed- 
ded in a linear space isomorphic to M'^' for some fc. 

Definition 1. Two generic histories H and H' , such that H{0) — H'{0), are said 
to be equivalent (in symbols H ~ H' ) when for every process Kp, with p > 0, 

a {Kp *H)^a{Kp* H') , z {Kp * H) ^ z {Kp x= iJ') , S {Kp *H)^S{Kp* H') . 
As a consequence, for p > 0, the condition of equivalence H ^ H' implies 

'■+°° / . . . X 

(GaW {s+p)W (s) + Gat. {s+p)v {s) + GaN {s + p) N {s)j ds ^ 

(g.w {s + p) W (s) + G„, {s + p) v' {s) + G^AT {s + p) N' (s)) ds 
[Gsw {s+p)W (s) -t- Gs. {s+p)iy (s) + GsN {s+p)N (s)) 
[Gsw {s + p) W {s) + Gsu {s + p) v' {s) + GsN {s + v) N' (s)) 







ds — 
ds 
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G,w {s +p)W (s) + G,^ {s+p)v (s) + G,N {s+p)N (s) j ds 

{g,w (s + p) W (s) + G,. (s + p) v' (s) + G.N {s + p) N' (s)) 
Let the distance d {H, H') be defined by 



ds. 



d{H,H') ~ sup 



+00 



G^w is + t)W (s) + G^^ is + t)i^ (s) + G^N {s + t)N (s) ) ds 
Gaw {s + 1) W' [s) + Ga^ [s + t) v' (s) + GaN (s + 1) N' (s)) ds + 
(^G,w {s + t)W (s) + G.^ {s + t),y (s) + G.n {s + t)N (s)) ds - 
Gzw {s + t) W (s) + G,^ {s + t) v' (s) + G.AT (s + t) N' (s)) ds + 



+ 00 



Gsw (s + t)W (s) + Gsu {s + t)v (s) + GsN {s + t)N {s) ds - 



- ^ (Gsw {s + t) W' (s) + Gsu {s + t) v' (s) + G5A' (s + t) N' (s)) ds | . 

It induces a pseudometric in the space F of histories and a metric in the quotient 
space T/ ^. By taking the hmit i — » in the expression above, one gets 

\a (H) - a {H')\ + \z {H) - z {H')\ + \S {H) - S {H')\ < d {H, H') + 
\{G„w (0) W (0) + G,, (0) u (0) + G^N (0) N (0)) - 
{Gaw (0) W (0) + Ga, (0) I/' (0) + GaN (0) A^' (0))| + 
\{G,w (0) (0) + G,, (0) 1/ (0) + G,N (0) iV (0)) - 
{G,w (0) W^' (0) + G,. (0) z/' (0) + G,N (0) TV' (0))| + 
I {Gsw (0) W (0) + G5. (0) V (0) + Gsw (0) (0)) - 
(Gsw (0) W^' (0) + Gsu (0) I/' (0) + GsN (0) iV' (0))| . 
Since two equivalent histories are characterized by identical initial values, namely 
W(0) = W (0), V (0) = t/' (0), N (0) = iV' (0), from previous inequality it follows 
that two equivalent histories determine the same macroscopic stress, microstress 
and suhstructural self-action. The proof of the analogous property for simple bodies 
in [21] is based on the use of a seminorm. 

Proposition 1. The pseudometric c?(-, ■) has the following properties: 

Contraction: for every p> Q 

d {Kp *H,Kp^ H') < d {H, H') . 

Fading memory: for every e > there exists r such that, for every p > r, 
one gets 

d {Kp *H,Kp* H') < e. 
Approachability: if Hp {p) = H' (0) , then 

lim d{H,Hp^H') =0, 

with Hp the process generated by H over [0,p). 
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Proof. Contraction arises directly from the definition. In fact, in writing explicitly 
d (Kp * H, Kp * H') by taking into account ()5.ip and the definition of the semimetric 
d (•, •), one manages integrals containing terms of the type Gab {s + t) B (s — r), 
with A and B running in {cr, z, S} and {W, respectively. Each of these terms 

is also equal to Gab [s + t + r) B (s) (see also dU). Consequently, by taking the 
superemum over p > one gets the contraction property straight away. 
Let now define 

M max < sup |W^'(s)|, sup |W^(s)|, sup |z^'(s)|, sup 11^(5)1, sup |A^'(s)|, sup |-/V(s)| 

Ls>0 s>0 s>0 s>0 s>0 s>0 

By taking into account that 

i-b 



Gab{s)B{s) ds 



< sup \B{s) 

s&[a,b] 



\Gab{s)\ ds, 



one obtains 



d {Kp *H,Kp* H') < 2M sup 

t>o 









GaW {s + t) 


Jp 





ds + 



+ 00 



+ 



G,w [s + 1) 
Gsw {s + 1) 



Ga^ [S + t) 

+ 00 



ds + 



ds 



Gsu {s + 1) 



GaN 


{s + t) 


ds+ 






r+oo 






ds -i 




GzN (s - 














r + co 






ds 4 


L 


GsN {s - 













ds- 



ds 



However, since it has been assumed that the maps s i — > \Gab \ (s) are integrable, 
there exists r such that, for r > m, one may find e > such that the right-hand 
side of the previous relation is lesser or equal to e {2AI) ^. Fading memory then 
follows. It also implies the property of approachability. In fact, from 

d{Hp * H', H) = d{Hp * H', Hp * W^) 

and fading memory, the approachability can be obtained by letting p to +00. □ 

Previous theorem suggests the following definition: 

Definition 2. H is said to be approachable from another history H' if there ex- 
ists a family of processes {p 1 — > Kp,p G K+) prolonging H' and such that Kp * H' 
converges to H' with respect to the pseudometric d{-, •) as p ^ +00. 

5.3. States and actions. The state space is identified here with the space of 
histories T endowed with the norm 



\H\\ 



\F\ 



In this way each state is defined to within an equivalent history. 



Definition 3. A function f : V 
f{H) = f{H'). 



is called a state function if H ^ H' implies 



Definition 4. A function a -.T x H ^ M. is called an action if 
(1): a is additive with respect to prolongations, namely 
a {K' x= K, H)^a {K' , K * H) + a {K, H) , 
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(2): the map a {K, •) : F — > M is continuous. 

Definition 5. The action a satisfies the dissipation property along H if, for 
every e > 0, there exists 5 > Q such that 

d{Kp* H,H) <S =^ a {Kp, H) > -e. 

Definition 6. Given a generic action a (•,•), a Junction / : F — > M is called a 
lower potential for a (■, •) if for every H and H' belonging to F and for every 
e > there exists 5 > Q such that 

f{H)-f{H')<a{K,H) + e, 

for every K such that d {K * H' , H) < S. 

The definition of state function comes from 21j - the difference is here only the 
extended meaning of the history - while the definitions of action and lower potential 
have been introduced in [12 . 

5.4. Work density in the bulk. At every x ^ B the work density w (H) is defined 

by 



w{H) 



a {H') ■ F{s) + z {H') ■ i^is) + S (H') ■ N{s) 



ds. 



Note the different algebraic sign in the analogous definition of w (H) given in [21] 
with reference to simple bodies, i.e. in absence of substructural interactions. 

Given an history H and a prolongation Kp — [Fp, Vp^ Np) of it, the work density 
over the prolongation Kp, indicated by w {Kp, H) is then defined by 

w {Kp, H) := w {Kp *H)- w{H) 

that is 

w {Kp, H) = r \a {k; * H) ■ Wp{s) + z {k; * H) ■ Vp{s) +S{k;*H)- Np{s) 



ds. 

An analogous power density over prolongations is defined in |21j but with reference 
to relative continuations (defined below). The use of relative prolongation implies 
the appearance of further terms in the explicit expression of w {Kd, H), terms due 
to the jump in H (0). 

Theorem 2. The work density over a prolongation is an action and for any fixed 
Kp the map w {Kp, ■) is a state function. 



Proof. To prove the latter property, first define 



M :— sup 



p 



Wp{s) ds, / \vp{s)\ ds, / Np{s) 



ds 



and remind that the manifold of substructural shapes is embedded in a linear space. 
It then follows that 

\w{Kp,H) -w{Kp,H')\ < Md{H,H') 

from which one realizes that w {Kp, •) is a state function. Previous inequality implies 
also that w {Kp, •) is also Lipschitz continuous. The additivity of w with respect to 
the processes is implied by the definition. □ 

The following lemma is a version of a proposition in [21] (see also [20]), the proof 
of which can be easily adapted to the present case. 



COMPLEX BODIES WITH MEMORY EFFECTS 



13 



Lemma 1. Assume that v can be freely selected in the linear space in which the 
manifold of suhstructural shapes is isometrically embedded. Given two different 
histories H — (W, v, N) and H' = {W ^ v\ N') and a process Lp (r) defined by 

Lp{r) := (^W{p) + -W'{0), ^i^{p) + -z^'(O), ^N{p) + -N'{0)) , 
\ P P P P P P J 

one gets 



lim d{Hp * Lp * H' . H) 

p— >+oo 



and 

w{Hp *Lp* H') = %u{H) + w{H') + 

+ Gab{^)B{^) ■ B{^) - Gab'{^)B'{0) ■ B'{0)). 

A,B A,B' 

Definition 7. A history H is said to be w— approachable from another history H' 
if H is approachable from H' and the sequence (p i — > Kp,p G M+) is such that the 
sequence p ' — > w (Kp, H) converges too. 

Theorem 3. // u can be freely selected in the linear space in which the mani- 
fold of substructural shapes is isometrically embedded, the space T is closed under 
w— approachability. 

The proofs of both the previous lemma and the approachabihty theorem foUow 
the same paths of the analogous results for simple bodies in [21] . The circumstance 
that V is now selected in a linear space implies just that, in re-following the path of 
the proofs in [5T], one needs only to consider the distance d (•, •) and the presence 
of the substructural terms. 

In order to prove the closure theorem (under w— approachability) the key point 
is the use of the last result of the Lemma. In fact, one replaces Hp * Lp with 

K2p ■■= {Hp - H' {0)1) * (Lp - H' (0)1) , 

where H' (0)^ is the constant history of value H' (0) and duration p, then one proves 
by Lemma that the work expended along the continuation K2p, namely w (i^2p, H) 
converges to w (H) plus the work done in the extreme retardation (see [H]) from 
H'{0) to iJ(oo). 

The work w helps also in characterizing the kernels in the constitutive expressions 
of the interaction measures. 

Definition 8. The relaxation functions s i — > Gab (s) are said to be dissipative if 

w{H)>0 

for any iJ G F such that F (oo) ^ 0, v (oo) = 0, iV (oo) ~ 0. 

Note that the requirement of the positivity of the work is weaker with respect 
5.5. ReleLxed work. 

Definition 9. For every pair of w— approachable histories H and H' , the relaxed 
work w^, (H) along H , starting from H' , is defined by 

Wh, (H) := mi I lim inf w{Kp, H') | /sTp G H, lim d{Kp *H',H) = q\ . 

I p^+oo p — > + oo J 
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Definition 10. If H is w— approachable from H' , w^r (H) is called the minimum 
work performed from H' to H while —w^' (H) is called the maximum recover- 
able work. 

The closure of the state space under approachabihty justifies the definition above. 
Moreover, an estimate follows: 

w'h' (H) <w{H) + ^{J2 Gab{^)B{^) ■ B(oo) - Gab'{^)B'{Q) ■ B'{Q)). 

AM A,B' 

It means that the relaxed power along H, starting from H' , is bounded from above 
by the power along H plus the difference of the powers under extreme retardation 
from H' (0) to H (oo). By restricting (•) to the histories K that are prolonga- 
tions of H itself, one gets an additional upper bound: 

w'h {K^ <w{k\h). 

Theorem 4. The following statements hold: 

(1) Both w^r (■) and w^.-j (H) are state functions. 

(2) (Sub-additivity.) For H, H' and H" histories such that w^i (H) > — oo and 
Whii {H') > — oo one gets the triangular inequality 

(3) (Lower semicontinuity.) // w^i (H) > — oo, w^i (•) is lower semicontinu- 

(4) (Dissipation inequality.) // w^, {H) > — oo and Kp is a process such that 

(Kp * H) > — oo, then 

w'h, {Kd * H) - w^H, (H) <w{Kp*H), 

moreover, if Kp is such that w^i {Kp * H) > — oo and w^^^^' (^) > 
then 

{H)<w{Kp,H'). 

The proof of the theorem above is essentially independent of the explicit expres- 
sion of the work. For this reason the proof of the analogous result in [2Tj applies 
providing one substitutes the norm used there with the distance defined above. 
Other results in [21] can be adapted here. Such results are listed below. Differ- 
ences rest essentially on (i) the presence of substructural terms, (ii) the use of the 
distance c? (•, •) and (iii) the use of strict continuations of histories, not the relative 
continuations used in |21| . The latter are indicated by R superposed to * and are 
defined by 

5m^.^.-/ Kp{s) + H{Q) 0<s<d 



Theorem 5. The following statements are equivalent: 

(1) w satisfies the dissipation property on all constant histories. 

(2) (H^) = along every constant history . 

(3) w^i (H) > for every history H. 



^Precisely, for any £ > there exists <5 > such that, given H, for any Hi such that d {Hi , H) < 
<5, one gets 

^H' (Hi)>w-^, {H)-e. 
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(4) For any history H one gets 

w{H) > Gab{oo)B{0) ■ B{0) - ^ Gab(oo)S(oo) • B{oo)). 

A,B A,B 

(5) The relaxation functions s i — > Gab (s) are dissipative. 

(6) w'^, {H) > for every pair of histories H and H' . 

(7) w satisfies the dissipation property for every history H. 

Theorem 6. The following statements hold: 

(1) For every history H the minimum power performed and the maximum re- 
coverable work from H to H are zero. 

(2) For every pair of histories H and H' 

(3) For every constant history H'^ , (•) and (-f^^) '^'^^ determined re- 
spectively by (•) and (O^); namely 

A,B 

w^H (H^) = w'h (0^) + ly2GAB {oo) BHi ■ BHf . 

^ A,B 

(4) The restriction of w"^ -^ (■) to constant histories is determined: for every pair 
of constant histories H\ and one gets 

^Hl (^2) = \{Y.GAB{<X>)B^yB^,^-Y,GAB{^)B^yB^,J. 

A,B A,B 

(5) For every pair of histories H and H' 



w^H' {H) > w^H' (h (0)t) = «;^,,(0t) + i ^ Gab{<^)B{Q) ■ B{0). 

A,B 

(6) For every H 



A,B 



(7) For every pair of histories H and H' 

^ A,B 

(8) For every pair of histories H and H' 

wl, {H) > w^H' (H) . 

(9) For every H 

w^,,(Ot) = inf w{K,H). 
Ken 

(10) For every constant history , the functional w^.^ (B^) upper semicon- 

tinuous. 

The symbols B^-i and Bu{q) used earlier indicate one of the entries of the list 
defining the state, evaluated along the constant history if^ or at H (0), respectively. 
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5.6. Free energies. The free energy can be defined in terms of actions (see fT2l). 
There are several possible free energies. Upper and lower bounds for their set can 
be determined. 

Definition 11. Every lower potential of w is called a free energy. 

The properties of the power discussed above allow one to prove the following 
theorem. 

Theorem 7. The following assertions hold: 

(1) Every free energy ip satisfies the dissipation inequality 

i}{K *H)-^{H) <w {K, H) 

for every H E T and every compatible K Cz H. Moreover, every l.s.c. 
function ip : T ^ R that satisfies the dissipation inequality is a free energy. 

(2) // the dissipation postulate is satisfied, then, for every H and H' belonging 
to T, the maps (•) and {H') are free energies. Moreover, for every 
free energy one gets 

-w^HiH')<^P{H)-i,iH')<w^H m 
for arbitrary histories H and H' . In particular, if there is H' and a family 
of free energies ip^ such that tp^ {H') = 0, the maps w^, (•) and """^[.^ {H') 
are the maximum and the minimum free energies in such a family. 

(3) Each free energy is a state function and, for every H , it satisfies the in- 
equality 

V'(i?(0)t) <^iH). 

In particular, the restriction of the free energy to constant histories is given 
by 

^ (i?t) _ ^ (ot) = IJ2Gab (oo) Bh, ■ Bh, . 

^ A,B 

The technique of the proof is strictly analogous (modulo the variations associated 
with the use of the metric d {■,■)) to the one used in [21] for an analogous result for 
simple viscoelastic bodies (see also [12]), so the details of the proof are not reported 
here. In |21j a weaker condition is adopted: discontinuity is admitted at under 
relative continuations. Here, the need of the use of the chain rule in an ensuing 
section suggests to avoid this discontinuity for the sake of simplicity. 

6. Characterization of the surface free energy in terms of the 

surface work 

The propositions presented so far can be extended in presence of structured 
discontinuity surfaces. In particular, it is assumed that across S the map v is 
continuous while the gradients W and N suffer bounded jumps. 

The attention is focused on constitutive relations of the type 



T(H) = Gtw(0)W(0) + GT^(0)i^(0) + Gtn(0)N(0)- 



+ 00 



G'tw(s)W(s) + Gt,{s)v{s) + Gtn(s)N(.s) 



ds, 
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S(e) = Gsw(0)W(0) + Gs.(0)i/(0) + Gsn(0)N(0)+ 



+ 



/• + 00 

Jo 



Gsw(s)W(s) + Gs4s)v{s) + Gsn(s)N(s) 



ds, 



3(H) = G3w(0)W(0) + G,^{OHO) + G,^{0)W{0) 



G3w(s)W(s) + Gi4s)u{s) + G3n(s)N(s)J ds, 
for the surface stresses T, S, and the surface self-action 3. 

Definition 12. Two surface histories M and M.', such that M (0) = H' (0) , are said 
to be equivalent if 

T(Kp*H) =T(Kp*H'), S(Kp*H) = S(Kp*H'), 3(Kp * H) = 3(Kp * H'), 
for every prolongation Kp. 

A distance (EI, EI') between surface histories can be defined by 



ds(H,H') := sup' 
t>o 



Gtw(s + t)W{s) - Gtw(s + i)W'(s) + Gj^s + t)u{s)- 



-Gt^s + t){sy{s) + Gtn(s + t)N{s) - Gtn(s + i)N'(s) 



ds 



+ 



+ 00 



G^w{s + i)W(s) - G3w(s + t)W'(s) + 63,(5 + t)u{s) 



-G3,(s + ty{s) + G3n(s + t)N{s) - G3n(s + t)N'{s) 



ds 



+ 



+ 



P + OQ 

Jo 



Gsw(s + t)W{s) - Gsw(s + t)W'(s) + Gsu{s + t)v{t)- 



-Gs^{s + ty{s) + Gm{s + t)n{s) - Gm{s + i)N'(s) ds } . 

It is a semimetric on the space of surface histories and a metric over the quotient 
space generated by the equivalence relation defined above. 

Proposition 2. The distance (■, ■) has the following properties: 

Contraction: for every r > 

dE(K^ *H',Kp *H) < rfE(H',]H[) Vr > 0. 

Fading memory: for e > there exists I such that, for every p > I, one 
gets 

(is(IKp *H',Kp *e) < e. 
Approachability: if Mp {p)~ = H' (0), then 

lim dE(Hp*e',H) = 0, 

p— >+oo 

with Mp the process generated by H over [0,p)- 
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The proof is analogous to the one of Proposition 1. In this case the constant M 
is the maximum of the suprema of the surface histories. 

Surface state functions and surface actions can be then defined. 

From Theorem 1 one reahzes that the surface work density is defined by 

:= J (^T • W + 3 ■(!>)+§• Nj dt + J {{a)m ■ [y] + {S)m ■ [v]) dt. 

It includes both peculiar surface interactions and traces of the bulk stresses at 
the discontinuity surface itself. Previous work on the instantaneous response of 
complex bodies with structured discontinuity surfaces j33| suggests that only a 
reduced surface work density is in strict connection with the surface energy 
density: 



:= f (t • W + 3 •(;>)+§• n) dt. 



The reduced surface work density over prolongations is defined by 

{Kp, H) := {Kp * H) - (H) , 

where Kp is the surface counterpart of Kp. 

By making use of the technique leading to Theorem 2, one may prove that 
is an action and the map H i — > {Kp, H) is a state function. 

A relaxed surface work can be then defined by 

(H) inf i liminf u)^ (K„, H') | K„ £ H, lim d(Kp * H', H) = 

As in the case of the bulk relaxed work, both (•) and wj^ (H) are state 
functions. Moreover, wj^ (•) is subadditive in the sense that, for histories H, H' 
and H" such that (H) > —oo and w^, (H') > — cxd, the inequality 

holds. If w"^ (H) > — cxD, (•) is lower semicontinuous. Under the same hypoth- 
esis, for every process Kp such that (Kp * H) > — cxo, the surface dissipation 
inequality 

{Kp * H) - (H) < (Kp * H) 

is verified. 

In summary, all the properties of w and hold also for and w^^ . The proofs 
can be constructed in the same way adopted in analyzing the power in the bulk. 

Definition 13. Every lower potential of is called a surface free energy. 

The theorem collecting the properties of the free energy in the bulk has its 
counterpart for the surface free energy. The proof is essentially the same. 

Theorem 8. The following assertions hold: 

(1) Every free energy (p satisfies the dissipation inequality 

(K =(= H) - (H) < WE (K, H) 

for every EI and every compatible K. Moreover, every l.s.c. function H i — > 
(p (H) G M that satisfies the dissipation inequality is a free energy. 
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(2) // the dissipation postulate is satisfied, then, for every pair of histories H 
and W, the maps Wjjp (■) and — (H') are free energies. Moreover, for 
every free energy (j) one gets 

-w^^ (HO < (H) - (H') < wlS" m 

for arbitrary histories H and H' . In particular, if there is H' and a family 
of free energies (p^ such that (p^ (H') = 0, the maps w'^ (•) and —w^^ (H') 
are the maximum and the minimum, free energies in such a family. 

(3) Every free energy is a state function and, for every H, it satisfies the in- 
equality 

^(h(0)^) < 0(HI). 

In particular, the restriction of the free energy to constant histories is given 
by 

ct> (Ht) - </> (Qt) = i ^ Gab {oo) Bnt • But • 

A,B 

A ranges in {T,3,S} and B in {W, i/,N}. 

7. The mechanical dissipation inequality and its consequences 

7.1. Mechanical dissipation inequality. In isothermal setting, the second law of 
thcrmodynamies in the form of Clausius-Duhem inequality reduces to a mechanical 
dissipation inequality. In Lagrangian representation, for any part b it reads 

±m{b,y,v)-Vr' (?)»<0. 

The functional \E' is the overall free energy of b along the fields y and v. it is the 
integral over b itself of the free energy density. If b is selected to cross E, so it 
is indicated by bs, both bulk and surface energy densities - the ones discussed 
previously - must be involved. Local forms of the mechanical dissipation inequality 
- local in the bulk and along E - can be obtained by exploiting the arbitrariness 
of the part considered. They are reported in summary here, written with reference 
to the infinitesimal deformation setting discussed in the earlier sections. The local 
form of the mechanical dissipation inequality in the bulk then reads 

i}-(j-W -z-u-S-N <Q, 

while the one at points over S is given by 

^ _ T • W - a-i/ - § • N < 0. 

7.2. Complex bodies with instantaneous elastic response. By borrowing 
terms from the mechanics of simple bodies, here bodies with instantaneous elastic 
response are the ones admitting constitutive structures of the type 

i^ = i;{H{t),H*), 

a=a{H{t),H'), 

z = z{H{t),H'), 
S = S{H{t),H'), 

in the bulk and 

= 0(H(t),H*), 
T = T(H(i),H*) , 
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3=3(H(<),rf), 
§ = S(H (t),H*) . 

The symbols H and H summarize the state. Precisely, H (t) — {W (t) , v (t) , N (t)) £ 
Max 3 X R'^ X Mkx3 is the state at the instant t while iJ* £ F is the past history of 
the state up to the instant t (in the notation used here H* is the graph of H from 
t to infinity, namely i/* (s) = H {t + s)). Analogous meaning can be attributed to 
H (t) and H*, namely H (<) = (W (t) , (t) , N (i)) £ Mgxs x R'^ x Mkx3 while H* is 
the past surface history and belongs to F. 

In the infinitesimal deformation setting treated here, it is assumed that the free 
energy density in the bulk is a quadratic form in the instantaneous values 
W (t), V [t) and N (t). It is also assumed that also the surface free energy 
density is a quadratic form in the instantaneous values W (i), v (<) and N (t). 

7.3. Chain rule. To exploit the local versions of the mechanical dissipation in- 
equality a chain rule must be used in evaluating the time derivative of the energy. 
Appropriate chain rules have been obtained in [37] and [2Qj, and can be adapted 
here. 

Consider a functional 

T : M3X3 X K'^ x A4x3 X F ^ M 

defined for every H {t) (or H(t)) in M3X3 x M'' x Mfcx3 and for every i/* (or 
H*) in F such that i?* (s) (or H* (s)) is in the open and connected subset U from 
M3X3 xR'' X Mfcx3, characterized by det{I + W) > (or det (/ + W) > 0), for 
almost s > (s is the time parametrizing the history 'prior' t - in the representation 
adopted here s > t). Assume that (i) T is continuously differentiable, (n) the 
function t 1 — > H [t) (or H(t)) with values in Li has two continuous derivatives 
t I — > H [t) and t \ — > H{t), and (m) for every t the past histories W and 
are in F. Under these conditions the function / (t) := T {H (t) , H*) (alternatively 
f (t) := (H (t) ,1131*)) is continuously differentiable and its time derivative is 

/ (t) = DT (H (t) ,H')-H (t) + ST [h (t) , iJ* 

(alternatively / (t) = DT (H {t) , H*)-]HI {t)+6J'{m (t) , H* |H*)), where DT {H {t) , H*) 
is a continuous functional taking values in T'^^^^IA for every fixed H (t) and iJ*, and 
6J- {H (t) , H*\K*) is a continuous scalar- valued functional depending linearly on 
if* and defined on the closed subspace of F spanned by the functions if* such 
that ii* (s) -I- if* (s) is in U for almost s > (analogous remarks hold also for 
DT (H (i) , H*) and 6T{m (t) , H*|IH*)). Fixed K\ an appropriate technical assump- 
tion is that ST {H (t) , H^\K*) is continuous in {H (t) ,H*). The proof of this chain 
rule can be found in pO] . 

7.4. Consequences of the mechanical dissipation inequality. As it is well 
known, to exploit the local version of the mechanical dissipation inequality one 
should have the possibility to select at will the instantaneous rate H (t) of the 
state. If this point is straightforward in the mechanics of (simple or complex) bodies 
without memory effects, some additional problems appear in presence of memory 
effects, due to the dependence of the interaction measures on the whole history of 
the state variables. The technique discussed in [20] to avoid these difficulties can 
be adapted here and is summarized in the following paragraphs, then it is applied 
to the case of complex bodies. 
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For instrumental reasons, it is useful to introduce a function / : M"'"^ R such 
that / (s) = for \s\ > 1, / (0) = 0, / (0) = 1. By using /, for any a G R+ and a 
fixed t S M+ one may define (see [Day], p. 91) varied histories 

[s] := H {s) + af (^^) (h - H {t)) , 

where W is a generic element from TH(t)U, namely W is a triple (V, v, T) of virtual 
rates of W , v and N . Essential properties of Ha (t) - in a sense the properties that 
suggest the definition of Ha {t) itself - are (i) Ha {t) = H (t), {ii) Ha {t) = 7i, {in) 
Ha {s) = H (s) for every s < t — a and s > t + a and (iv) for a sufficiently small 
Ha (•) meets the hypotheses of the chain rule and both Ha (•) and Ha (•) converge 
in norm respectively to H (•) and H (•) as q ^ 0. 

A similar definition can be adopted for the surface history (s) so that one 
gets 

H„ (s) := m (s) + af (^^^ (sj-U (i)) , 

where, now, ij is a generic element from T]Hi(t)Z//, namely Sj is a triple (QJ, v, &) of 
virtual rates of W, i' and N. 

By making use of the chain rule and substituting H (s) with Ha (s), from the 
local mechanical dissipation inequality in the bulk one gets 

{dwit)i^ {H {t) ,H*)-a [H {t) ,H*))-V + {H (t) ,H')-z {H {t) ,H'))- v+ 

+{dNit)^ {H (i) ,H')-S {H it) ,H')).r + Si,(H (t) , H'\H') < 

as a — > 0, an inequality holding for all choices of the triple (V, v, T). The arbitrari- 
ness of (V, V, T) implies that in the bulk 

a{H{t),H')=dwit)^{H{t),H'), 

z{H{t),H*)=d,^t)i:{H{t),H*), 
S{H{t),H')^dN(^t)4'{H{t),H'), 

Si; (^H{t),H'\H'^ < 0. 
An analogous result hold along the surface S where, locally one gets 
(5w(t)<^ (H (t) , H*) - T (H (t) , H*)) • 5J + (H (t) , H*) - 3 (M (t) , H*)) • v+ 

(H (t) , H*) - 5 (M (t) ,m*)) ■ & + S(l) (m (t) , H*|IHI*) < 

as a — ^ 0, an inequality holding for all choices of the triple (2J, v, <3). The arbitrari- 
ness of (2J, u, 6) implies that, along the surface, 

T(H(t),H*) =%(t)(/.(H(f),H*), 
3(H(i),H*) =5,(t)0(H(f),H*), 

5(H(t),H*) = aN(t)(A(iHi(t),e*), 

5(j!)(^H(i),e*|lHI*) < 0. 
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